1902.] 


On Spherical Harmonic Analysis. 


97 


“ On some Definite Integrals and a Hew Method of reducing a 
Function of Spherical Co-ordinates to a Series of Spherical 
Harmonics/’ By Arthur Schuster, F.R.S. Received May 
30,—Read June 5, 1902. 

(Abstract.) 

The expansion of a function f(0) of an angle 0 varying between 0 
and 7r in terms of a series proceeding by the sines of the multiples of 
0 depends on the fundamental theorem, 


| sin pO sin qO clQ = 0, 
o 

where p and q are integer numbers not equal to each other. Simi¬ 
larly if P n denotes the zonal harmonic of degree n 9 g = cos 0 , and 

sl<rT> 

Q* = sin 0, 0 , 

dpn 

the expansion of a function of 0 in terms of a series of the functions 
Q£ depends on the corresponding theorem, 

+i 

| QSQ?^ = o, 

~i 

where i and n are two integer numbers not equal to each other. In 
many practical applications a continuous function is given by means 
of its numerical values at certain points, e.g ., for equidistant values 
of 0. 

Such cases present no difficulty when Fourier’s analysis is to be 
employed, because there is in that case a summation theorem exactly 
corresponding to the above integration theorem. If 0 be replaced by 
pir/n, where p takes successively the values 1 , 2, 3 ... , the equa¬ 
tion 

p = n —1 

2 sin (ppir/n) sin (pq^/n) — 0 
o = o 

will hold true. This allows us to determine the coefficients in the 
case of problems in which discontinuous values of the function at 
equidistant points are known (e.g., hourly readings of temperature or 
barometric pressure). If we assume that all Fourier coefficients beyond 
the nth vanish, n equations are obtained, each of which only contains 
one of the unknown quantities. 

If it is desired to expand a function in terms of cosines, a slight 
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modification must be introduced, the summation theorem in that case 
being 

p — n—l 

b + b cos pir COS (J7r -f- 2 COS (ppir/n) COS (pcprhi) = 0, 

p = i. 

the first and second terms representing half the value of the product 
for p == 0 and p — n respectively. 

There is no corresponding summation theorem in the case of the 
functions Qj, and the application of the method of least squares leads 
to a series of normal equations, each of which contains all the other 
coefficients. This has been one of the great practical difficulties in 
obtaining an expression for the series of spherical harmonics for the 
earth’s magnetic potential. 

F. E. Neumann has tried to overcome the difficulty by calculating 
coefficients cq, . . . a q in such a way that 

V 0pQJ fap)QfO*p) » 0, 

p = i 

Here /xi, g 2 * • • gg are the quantities for which the values of the 
function to be represented are known. Neumann’s process is equiva¬ 
lent to attaching weights proportional to a p to the different obser¬ 
vations, a proceeding against which theoretical objections might be 
urged. 

2. The expansion in terms of a series of cosines and sines being so 
much easier than the direct expansion in terms of a series of the 
functions Q£, I have endeavoured to obtain the latter series by means 
of the former. 

It is well known that a function of an angle 6 , which is confined 
to the values lying between 0 and t r, may be put either into the form 

a 0 + tii cos 0 + a.o cos 26 + . . . a p cos pO +. . ., 
or into the form 

b'i sin O + h, sin 20 -f- . . . b 2 j sin r p0 + . . . 

The reduction to the series of spherical harmonics is accomplished 
by calculating and tabulating the coefficients in the series 

cos 1*9 = :A£<£ +a; +1 Q' +1 + . . . azq: 

sin P 6 =B'QJ+B' +1 QJ +1 + . . . B'QJ + . . . 

The choice between the cosine and sine series is open to us, but it 
appears that great simplicity is gained by taking the former series 
when <r is odd and the latter when <r is even. For in that case the 
coefficients A,f and will all vanish, as long as n is smaller than p. 
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When it is therefore desired to retain only terms as far as the nth 
degree, the Fourier coefficients need only be calculated as far as 
p = n+l. The position of the earth’s magnetic axis, e.g ., only 
depending on the terms of the first degree, is completely determined 
by the coefficients b 2 for cr = 0 and a 2 for o- = 1. 

3. The symbolical representation of the results of this paper is much 
facilitated by the introduction of a separate symbol for the product of 
alternate f actors, n . n - 2 . n - 4 ... 1, if % be odd, or n.n- 2 
. . . 2 if n be odd. I propose to write n !! for such products, and if a 
name be required for the product to call it the “ alternate factorial ” 
or the “ double factorial.” Full advantage of the new symbol is only 
gained by extending its meaning to negative values of n. Its complete 
definition may then be included in the equations 


n !! = n(n- 2 )!!, 1!! - 1, 2 !! = 2 . 

From this we may derive when n is negative and odd 


<-l ) 7 


1 

- n 2 ) \i 5 


while for n negative and even, the factorial becomes infinitely large. 

4. The calculation of the factors Aj and Bj depends on the values 
of the definite integrals 

-i +i 

| Qj cos pOclfi, j Qj sin pOdp , 

-i -i 


and these may be made to depend on the values of the integrals 



It is proved that 


+i 

and j /xQj sin x 0dy . 

-i 


J Q» sin '= c /,, vz-- T-'- n! !■; v:r; ' v w rn . if n - be even, 


-1 
+ 1 


/xsin x 6dy = c 


(n - cr) !! (a* — A — 2)!! (71+A-fl)!! 

0 if n - o- be odd, 

(n + <r) ![(cr-h A)!! (n - A - 3) 1! 


(n- 


1)!! (cr — A - 2)!! (n+ A + 2)!! 


0 


, if n - a- be odd, 
if n - cr be even. 


The factor c is equal to 2 or to 7r according as cr + A is even or odd. 

5. The integrals 



100 


Prof. A. Schuster. [May 30, 


are obtained in the form of a series having a finite number of 
terms. 

6. To find 

+i +i 

j" sin pOdg and [ QJ cos p9d[i, 

-i ~i 


we may either express Qj or the trigonometrical functions in terms of 
a series of powers of sin 0 . The second alternative leads to results 
which in general are more convenient. 

If we put 


Co - 1; 


Ci =■ p; 

C A 


Vi = p • I; 
£ 


c 3 


'•P 


p -1. p+ i 

■ l 72 73“" ; 


p 0?+x~2)!! m 

X ! (p - X) 1! J 


B 0 = 1; Bi = p; B 3 


p - 1 .p + 1 # 

T72 ; 




__ p-2 .p .p + 2 

^ 1 ,2.3“ 


1 (p + X-1)!! 

X ! (j? - X ~ 1) !! ’ 


we find if <r be even, p odd, and n even, 

j*Q sin pOdji — w - <fCi.cr-l.cr+l 

J •» 1 r (n - a-)!! (w + 2) !! I 

~ <r — 3.(r — 1 .ct- + 1.ct- + 3 ^ cr — 5. <r — 3.(r — l.<r+l.o*+3.<r + 5 

— L/3------+ O5 - 


n - 3. ?i + 4 


n- 3.n-5.n + 4c.n + § 


and if o- be even, p even, and n even, 


(n + (r)!! (n — 4)!! 
(n — cr — I) !! (pi + 3) 


JJLf 

3)!! I 


Bi.<r~l.<r+1 — B 3 


O' — 3 . (T — 1 . (T + 1 . O* — 3 

n — 4 .11 + 5 


+ B f 


cr — 5.cr—3.cr — l.cr+l.cr+3.cr+5 
n - 4. n - 6. n + 5. n + 7 




If n + cr -f-_p be even, the integral is zero. 
Similar equations are obtained for 


+ 1 

| Qj cos pQd[ a. 
-1 


7. The final results are expressed as follows :— 

If gZ denote the coefficient of cos ar<j> in the series of spherical 
harmonics when R 0 " = ££Qj, and t% is a numerical coefficient, </> being 
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the longitude, the results of the investigation may he put into the 
form 

P — ?t+l 

(f n — 2 (f r *, when n is even, <r odd, and p odd. 

p — o 

P = 91+1 

= 2 < K » n » odd, ,r » » P even. 

p=i 

p = 91 + 1 

= 2 „ n „. even, or even, „ p odd. 

v -1 

p — w + l 

= 2 % *» ». n ». odd, <r „ ,, p even. 

:p =2 

In these equations the factors b% are the coefficients of the 
Fourier series (see § 2), and the quantities rj, mj, sj are 
numerical quantities, which (as well as their logarithms) are gisren in 
tables at the end of the paper as far as n = 12, o- = 12, p — 12. By 
means of these tables the numerical work is reduced to a minimum, and 
the coefficients of the series may be obtained as far as terms of the 
12th degree. 

8. The proposed method is specially adapted to deal with problems 
like that of terrestrial magnetism, in which the function to be 
obtained as a series of spherical harmonics is not given directly, 
but by means of its differential coefficients. The force directed to 
the geographical north may by Fourier’s analysis be obtained as 
a sum, the terms of which have the form cos o-<p cos pO , and sin <rcf> 
cos p9 when cr is even, and the form cos cr<£ sin ^90, sin o-cj> sin pO wheno- 
is odd. Integrating with respect to (9, the magnetic potential is obtained 
in a form such that the transformation into the series of spherical 
harmonics may be proceeded with. A separate expression of the 
magnetic potential is derived from the force directed to the geographical 
east. 



